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The massive non-Abelian gauge fields are quantized Lorentz-covariantly 
in the Hamiltonian path-integral formalism. In the quantization, the Lorentz 
condition, as a necessary constraint, is introduced initially and incorporated 
into the massive Yang-Mills Lagrangian by the Lagrange multiplier method so 
as to make each temporal component of a vector potential to have a canon- 
ically conjugate counterpart. The result of this quantization is confirmed by 
(^-^ ■ the quantization performed in the Lagrangian path-integral formalism by ap- 

plying the Lagrange multiplier method which is shown to be equivalent to the 
, Faddeev-Popov approach. 
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I. INTRODUCTION 



It was shown in [1,2] that the massive non-Abelian gauge fields can well be quantized 
in the Hamiltonian path-integral formalism following the procedure initiated in [3,4] for 
the massless gauge fields. In this quantization, the massive Yang-Mills Lagrangian density 
written below was chosen to be the starting point. 

C = ~F a ^F^ + lm 2 A a ^ (1.1) 

. where are the vector potentials for a massive gauge field, 

O : F; u = d u A a u - d v A« + g! ahc A\Al (1.2) 



are the field strengths and m is the mass of gauge bosons. Working with the Lagrangian in 
equation (1.1), the massive non-Abelian gauge fields appear to be second class systems. In 
the quantization, it is necessary to introduce a primary constraint of second class defined by 

n S^) = ^ = ° (1-3) 



dA% 



1_h ' which was incorporated into the Lagrangian (1.1) by the Lagrange multiplier method. The 
result of the quantization amounts to the one given in so-called unitary gauge. An alternative 
quantization within the Hamiltonian path-integral formalism was subsequently performed in 
[5] by the approach proposed initially in [6] . In this approach, the second class non-Abelian 
system is converted to a first class system by enlarging the phase space with introducing a 
series of extra fields. The results of this quantization can be compared with those obtained 
early in the Lagrangian path-integral formalism by using the gauge-invariant Stiickelberg 
Lagrangian [7-11]. Due to the introduction of the extra field variables which have a corre- 
spondence with the Stiickelberg scalar functions, the quantized result looks much compli- 
cated. 
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In this paper, a Lorentz-covariant quantization of the massive non-Abelian gauge fields will 
be carried out within the framework of Hamiltonian path-integral along the line described 
in our previous paper for the quantization of massless non-Abelian gauge fields [12]. The 
essential point of this quantization which is different from the previous is that the Lorentz 
condition 

tp a = 0M£ = (1.4) 

, as a necessary constraint, is introduced from the beginning and imposed on the massive 
Yang-Mills Lagrangian. This treatment is consistent with the fact that a massive gauge 
field has only three polarization states and can completely be described by the Lorentz- 
covariant (four-dimensionally) transverse part of the vector potential, A^(x). Whereas, the 
Lorentz-covariant longitudinal part of the vector potential, A"^ , appears to be a redundant 
unphysical variable which must be constrained by introducing the Lorentz condition which 
implies A a ^ — 0. Conventionally, the Lorentz condition is viewed as a consequence of the 
field equations of motion [13] 

d»F» v + m 2 Al=j a v (1.5) 

where j£ is the current generated by the gauge field itself. The argument of this viewpoint 
is as follows. When we take divergence of the both sides of equation (1.5) and notice the 
current conservation, it is found that 

m 2 0"A° - (1.6) 

Since m ^ 0, the above equation leads to the Lorentz condition. This seems to imply that 
the Lorentz condition has already been included in the massive Yang-Mills Lagrangian. If so, 
when the Lagrangian is written in the first order form, we should see a term in the Lagrangian 
which is given by incorporating the Lorentz condition by the Lagrange multiplier method. 
Nevertheless, as will be seen in the next section, there is no such a term to appear in the 
Lagrangian. Therefore, the viewpoint stated above is not reasonable. The correct procedure 
is to treat the Lorentz condition as a primary constraint imposed on the massive Yang-Mills 
Lagrangian. In this case, due to the Lorentz condition introduced, equation (1.6), as a 
trivial identity, naturally holds, exhibiting the self-consistency of the theory. This procedure 
coincides with the aforementioned procedure of the quantization performed in [1,2,5] where 
the condition in equation (1.2), as a primary constraint, is necessarily introduced and is also 
consistent with the conventional canonical quantization. In the latter quantization, to derive 
the free propagator of the gauge boson, one only uses the transverse free field operator [13] 

A° fl (x) = / Ve^(fc)[a c A (fc)e-^+a- + (fc)e fe ] (1.7) 

J (27r)VM*0 a=i 

where (fc) are the polarization vectors satisfying the transversality condition 

fc"e£(fc)=0 (1.8) 

which follows directly from the Lorentz condition. 

We have a great interest to note that in the physical subspace defined by the Lorentz 
condition, i.e., spanned by the transverse vector potential A^i 1 , the dynamics of massive 
gauge fields is gauge-invariant in contrast to that the massive Yang-Mills Lagrangian itself 
is not gauge-invariant. In fact, when we make an infinitesimal gauge transformation[13 ]: 
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SA; = Df9 b (1.9) 

where 

Df = 5 a %-gf abc A* (1.10) 

to the action given by the Lagrangian in equation (1.1), noticing the identity f abc A a t t A b fi = 
and applying the Lorentz condition, it can be found that 

SS = -m 2 J d i x6 a d^Al = (1.11) 

Alternatively, the gauge-invariance may also be formulated by means of the Lagrangian 
represented in terms of the transverse fields 

C = -\F a T ^F^ v + \m 2 A a T »A^ (1.12) 

which is obtained from equation (1.1) by applying the solution of equation (1.4): A a £ = 0. 
Under the gauge transformation taking place in the physical subspace 

SA^ = D^9 b (1.13) 

where Dj? is defined as that in equation (1.10) with A c ^ being replaced by Aj. , it is easy 
to prove that the action given by the Lagrangian (1.12) is invariant 

5S = -m 2 J d A xQ a d ti M Fll = (1.14) 

where the transversality condition: d^Aj, — 0, as an identity, has been noticed. The gauge- 
invariance of the action in the physical subspace suggests that the quantum non-Abelian 
gauge field theory may also be set up on the basis of gauge-invariance principle as will be 
specified in section 4. 

The remainder of this paper is arranged as follows. In section 2, we will start from the 
Lagrangian written in the first order formulation. This Lagrangian is given by incorporating 
the Lorentz condition into the massive Yang-Mills Lagrangian by the Lagrange undetermined 
multiplier method and will be used to derive canonical equations of motion for the fields. It 
will be demonstrated that the equations of motion derived are self-consistent and complete 
for describing the field dynamics. In scction3, the quantization of the massive gauge field 
will be Lorentz-covariantly performed in the Hamiltonian path-integral formalism. The re- 
sult of this quantization is different from those given in [1,2,5]. In section 4, for comparison, 
the quantization will be alternatively carried out in the Lagrangian path-integral formalism 
by applying the Lagrangian multiplier method and the gauge-invariance principle. It will 
be shown that the Lagrangian multiplier method is equivalent to the Faddccv-Popov ap- 
proach [14]. The last section serves to make conclusions and comments on the problems of 
renormalizability and unitarity of the theory. 



II. FIRST ORDER FORMULATION AND EQUATIONS OF MOTION 

According to the general procedure, the Lorentz condition (1.4) may be incorporated into 
the Lagrangian (1.1) by the Lagrange undetermined multiplier method to give a generalized 
Lagrangian. In the first order formalism, this Lagrangian is written as [3,4] 
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C = \F a ^F; v - ±F a ^(d„A a u - d v A^ + gf abc AlAl) + \ m 2 A a » A^ + A a 5^ (2.1) 

where A^ and F^ v are now treated as the mutually independent variables and A a are chosen 
to represent the Lagrange multipliers. Using the canonically conjugate variables defined by 

TT M - ^ _ F a i \or r F ko = E k> lf A 1 = k = i; 2 ' 3 ; (o o) 

H \ ) - fj.0 ' <Vo - i X a = _ E a^ j f [L = Q l^J 

the Lagrangian (2.1) may be rewritten in the canonical form 

C = ET^A" + A a a C a - E^ a - H (2.3) 

where 

C a = dt i E a + gfabc^ck + ^ ^ 

n = \{EZf + \{F^ + \m 2 [{Atf + (At) 2 ] (2.5) 

here = (Eq,E%) is a Lorentz vector, H is the Hamiltonian density in which F°- are 
defined by 

F* = diA" - djAI + gr hc A\A) (2.6) 

In the above, the four-dimensional and the spatial indices are respectively denoted by the 
Greek and Latin letters. From the stationary condition of the action constructed by the 
Lagrangian (2.3), one may derive the equations of motion as follows 

A% = d k A« + gf abc A\Al - El (2.7) 
K = d l F? k + gf abc (E b A c Q + F b t A") + m 2 A% + d k E« (2.8) 

C a (x) = d^El + gf abc A b k E ck + m 2 Al = (2.9) 

and equation (1.4) where k — 1,2,3. Equations (2.7) and (2.8) act as the equations of 
motion satisfied by the independent canonical variables At and E^ (k = 1, 2.3) which precisely 
describe the three degrees of freedom of polarization for a massive gauge field with a given 
group index, while, equations (1.4) and (2.9) can only be regarded as the constraint equations 
obeyed by the constrained variables Aq and Eq because in these equations, there are no time- 
derivatives of the dynamical variables Af, and E^. Equation (2.3) clearly shows that these 
constraints have already been incorporated in the Lagrangian by the Lagrange undetermined 
multiplier method. Especially, the Lagrange multipliers are just the constrained variables 
themselves in this case. It is clear to see that in equations (2.7)-(2.9) and (1.4) there are 
altogether eight equations for a given group index. They are sufficient to determine the eight 
variables including three pairs of the dynamical canonical variables At and E^(k = 1,2.3) 
and one pair constrained variables Aq and Eq , showing the completeness of the equations. 

Along the general line by Dirac [3] , we shall examine the evolution of the constraints (p a 
and C° with time. Taking the derivative of the both equations (1.4) and (2.9) with respect 
to time and making use of the equations of motion : 
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which are obtained from the stationary condition of the action given by the Lagrangian (2.3), 
one may derive the following consistency equations [2,3] 



{ff, ¥>"(*)} + J d 4 y[{^(x),C b (y)}A b (y) - {<p°(x), V b {y)}E b (y)\ = 



(2.12) 



{H, C a (x)} + J d 4 y[{C a (x), C b (y)}A b (y) - {C a (x), ^ b (y)}E b (y)} = (2.13) 

where equations (1.4) and (2.9) have been used. In the above, H is the Hamiltonian defined 
by an integral of the Hamiltonian density shown in equation (2.5) over the coordinate x and 
{F, G} represents the Poisson bracket which is Lorentz-covariantly defined as 

Id 4 — 5G — SG 

The Poisson brackets in equations (2.12) and (2.13) are easily calculated. The results are 



{C a {x),y b { y )}=Df{x)d»8\x-y) 


(2.15) 


W a (x),<p»(y)} = 


(2.16) 


{C a (x),C b (y)} = m 2 [gf abc A c (x) - 25 ab dS]5 4 (x - y) 


(2.17) 


{H,<p a (x)} = d k x E a k {x) 


(2.18) 


{H, C a (x)} = m 2 [dSA a (x) + d%A%{x)] 


(2.19) 



where D® b (x) was defined in equation (1.10). It is pointed out that by the requirement of 
Lorentz-covariance, in the computation of the above brackets, the second term in equation 
(2.9) has been written in a Lorentz-covariant form gf abc A bf *E^. We are allowed to do it 
because the added term gf abc AQEQ only gives a vanishing contribution to the term AqC 11 in 
equation (2.3) due to the identity J^AqAq = 0. Particularly, the determinant of the matrix 
which is constructed by the Poisson bracket denoted in equation (2.15) is not singular. This 
indicates that equations (2.12) and (2.13) are solvable to determine the Lagrange multipliers 
Aq(x) and Eq(x). There is no necessity of taking other subsidiary constraint conditions 
into account further. This reveals the consistency of the constraints (1.4) and (2.9). On 
substituting equations (2.15)-(2.19) into equations (2.12) and (2.13), we find 

n x A a (x) gr bc d^A b (x)A^x)} - d k x E a k {x) = (2.20) 
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[5 ab n x -gf<*°A°(x)d£]E%(x) = 



(2.21) 



These equations are compatible with the equations (1.4) and (2.7)-(2.9). In fact, as can 
easily be verified, equations. (2.20) and (2.21) may directly be derived from equations (1.4) 
and (2.7-(2.9). In addition, we mention that if the Hamiltonian density is defined by H = 
H — A^C a + E^ifi 11 , the equations (2.20) and (2.21) can also be obtained from the equations 
{H, ip a (x)} = and {H, C a (x)} — respectively where H is the Hamiltonian defined by an 
integral of the H over the coordinate x. 

III. QUANTIZATION IN THE HAMILTONIAN PATH-INTEGRAL FORMALISM 

This section serves to formulate the quantization performed in the Hamiltonian path- 
integral formalism for the massive non-Abelian gauge fields. In accordance with the general 
procedure of the quantization [4] , we firstly write the generating functional of Green's func- 
tions via the independent canonical variables which are now chosen to be the transverse 
parts of the vectors A^ and E° 

Z[J] = 1 J D(A a /, E^)exp{i J d 4 x[E a /A a T ^ - H*(A a /, E a /) + J^A^]} (3.1) 

where 7i*(Ajf , E^t 1 ) is the Hamiltonian which is obtained from the Hamiltonian (2.5) by 
replacing the constrained variables A a ^ and E^ 1 with the solutions of equations (1.4) and 
(2.9). As mentioned before, equation (1.4) leads to A a ^ = 0. Noticing this solution and the 
decomposition E a ^(x) = E^ + E^(x), when setting E^(x) = d£Q a (x) where Q a (x) is a 
scalar function, one may get from equation (2.9) an equation obeyed by the scalar function 
Q a (x) for a given group index 

K ab (x)Q\x) = W a (x) (3.2) 

where 

K ab {x) = 6 ab n x - gf abc A^{x)d^ (3.3) 

and 

W a (x) = gr bc E b /(x)A c Tt ,(x) - m 2 Af(x) (3.4) 

With the aid of the Green's function G ab {x—y) (the ghost particle propagator) which satisfies 
the following equation 

K ac (x)G cb (x -y) = S ab S 4 (x - y) (3.5) 
one may find the solution to the equation (3.2) as follows 

Q a (x) = J d 4 yG ab (x - y)W b {y) (3.6) 

With the expressions given in equations (3.4) and (3.6), we see that the E'^ 1 (x) is a com- 
plicated functional of the variables A^f and so that the Hamiltonian H* (A^t 1 , E^) is 
of much more complicated functional structure which is not convenient for constructing the 
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diagram technique in perturbation theory. Therefore, it is better to express the generating 
functional in equation (3.1) in terms of the variables A® and E^. For this purpose, it is 
necessary to insert the following delta-functional into equation (3.1) 

5[Af]8[E a L ^ - E^(A^,E^)] = detM 5[C a ]5[^ a ] (3.7) 

where M is the matrix whose elements are 

M ab {x,y) = {C a {x),y b {y)} (3.8) 

which was given in equation (2.15). The relation in equation (3.7) is easily derived from equa- 
tions (1.4) and (2.9) by applying the property of delta- functional [4]. Upon inserting equation 
(3.7) into equation (3.1) and utilizing the Fourier representation of the delta-functional 

S[C a ] = J D( V a /2iryJ d4xvaca (3.9) 

we have 

Z[J] = iy D(A^EZ,ri a )detM6[d»A«}exp{i J d 4 x[E a »A« 

+r] a C a -H(A a ^E%) + J^A^]} (3.10) 

In the above exponential, there is a Eft -related term Eft {doA% — dorj a ) which permits us to 
perform the integration over Eft, giving a delta-functional <5[<9oAg — doi] a ] = det\do\~ 1 S[AQ — 
rj a }. The determinant det | c?o | 1 , as a constant, may be put in the normalization constant 
N and the delta- functional S[Aq — rj a ] will disappear when the integration over r/ a is carried 
out. The integral over E% is of Gaussian-type and hence easily calculated. After these 
manipulations, we arrive at 

Z[J] = 1 y D(A%)detM5[d*A$exp{i j d^-F^F^ 

+ l m *A a i>A° + J^Affl (3.11) 

When employing the familiar expression [4] 

detM = J D{C a ,C a )e i f dixdiyea{x)Mab{x ' v)cb{y) (3.12) 

where C a (x) and C a (x) are the mutually conjugate ghost field variables and the following 
limit for the Fresnel functional 

5[d»A a u ] = lim C[a}e-^f d4x(d " A ^ 2 (3.13) 

where C[a] <~ Ilajl^^) 1 ^ 2 an< ^ supplementing the external source terms for the ghost fields, 
the generating functional in equation (3.11) is finally given in the form 

Z[J,l£] = ^ j D(A^,C a ,C a )exp{i J d 4 x[C eff + J^Al+rC a + C a C}} (3.14) 
where 
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C eff = ~F a ^ v F^ + ^m 2 A a ^A« - ^(^A^f ~ d tl C a D ab C b (3.15) 

which is the effective Lagrangian for the quantized massive gauge field in which the first two 
terms are the Yang-Mills Lagrangian , the third and fourth terms are the so-called gauge- 
fixing term and the ghost term respectively. In equation (3.14), the limit a — > is implied. 
Certainly, the theory may be given in general gauges (a / 0). In this will be seen 

in the next section, the ghost particle will acquire a spurious mass /j, = ^am. 

IV. QUANTIZATION IN THE LAGRANGIAN PATH-INTEGRAL FORMALISM 

To help understanding of the result of the quantization given in the preceding section, in 
this section, we attempt to quantize the massive non-Abelian gauge fields in the Lagrangian 
path-integral formalism. For later convenience, the massive Yang-Mills Lagrangian in equa- 
tion (1.1) and the Lorentz constraint condition in equation (1.4) are respectively generalized 
to the following forms 

C X = -\F a ^F; u + \m 2 A a »Al - \a{X a f (4.1) 

and 

0M£ + aX a = (4.2) 

where X a (x) are the extra functions which will be identified with the Lagrange multipliers 
and a is an arbitrary constant playing the role of gauge parameter. Now, according to the 
general procedure for constrained systems, equation (4.2) may be incorporated into equation 
(4.1) by the Lagrange multiplier method, giving a generalized Lagrangian 

C x = ~F a » v F^ + ^m 2 A af *A a ^ + A°0"i4° + ^a(X a ) 2 (4.3) 

This Lagrangian is obviously not gauge-invariant. However, for building up a correct gauge 
field theory, it is necessary to require the dynamics of the gauge field to be gauge- invariant. 
In other words, the action given by the Lagrangian (4.3) is required to be invariant under the 
gauge transformations shown in equation (1.9). By this requirement, noticing the identity 
f abc A afJ -A b l = and applying the constraint condition (4.2), we have 

<SS A =~f d A xd v A a v {x)di i {Vf{x)d b {x)) = (4.4) 

where 

V?(x)=5° b £-&Z + D?(x) (4.5) 

in which /j, 2 = am 2 and D® b (x) was defined in equation (1.10). From equation (4.2) we see 
\d v A1 = — X a 7^ 0. Therefore, to ensure the action to be gauge-invariant, the following 
constraint condition on the gauge group is necessary to be required 

dH(Vf(x)e b (x))=0 (4.6) 
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These are the coupled equations satisfied by the parametric functions 9 a (x) of the gauge 
group. Since the Jacobian is not singular 



detM ^ (4.7) 

where 



M (x, y) = 



S6 b (y) 

= 5 ab {u x + v?)6\x - y) - gf abc dH(Ae(x)5*(x - y)) (4.8) 

the above equations are solvable and would give a set of solutions which express the functions 
6 a (x) as functionals of the vector potentials A^(x). The constraint conditions in equation 
(4.6) may also be inserted into the Lagrangian (4.3) by the Lagrange undetermined multiplier 
method. In doing this, it is convenient, as usually done, to introduce ghost field variables 
C a {x) in such a fashion [13] 

e a ( x ) = C,C a {x) (4.9) 

where £ is an infinitesimal Grassmann's number. In accordance with equation (4.9), the 
constraint condition (4.6) can be rewritten as 

d»(VfC b ) = (4.10) 

where the number £ has been dropped. This constraint condition usually is called ghost equa- 
tion. When the condition (4.10) is incorporated into the Lagrangian (4.3) by the Lagrange 
multiplier method, we obtain a more generalized Lagrangian as follows 

C x = ~F a ' a, F^ v + ^m 2 A a ^Al + A o 0"A° + ^a(X a ) 2 + C a d^{V ab C b ) (4.11) 

where C a (x), acting as Lagrange undetermined multipliers, arc the new scalar variables 
conjugate to the ghost variables C a (x). 

At present, we are ready to formulate the quantization of the massive gauge field . As we 
learn from the Lagrange undetermined multiplier method, the dynamical and constrained 
variables as well as the Lagrange multipliers in the Lagrangian (4.11) can all be treated as 
free ones, varying arbitrarily. Therefore, we are allowed to use this kind of Lagrangian to 
construct the generating functional of Green's functions 

Z[J a ",?,?] = ± J D(A«,C a ,C a ,\ a )exp{i J d 4 x[C x (x) + J^(x)A^(x) 

+T(x)C a (x) + C a (x)C(x)}} (4.12) 

where D(A" , • • • , A a ) denotes the functional integration measure, J®,K a and K a are the 
external sources coupled to the gauge and ghost fields and N is a normalization constant. 
Looking at the expression of the Lagrangian (4.11), we see, the integral over X a (x) is of 
Gaussian-type. Upon completing the calculation of this integral, we finally arrive at 

Z[^,f,^] = l J D(A^,C a ,C a ,)exp{i J d 4 x[£ e// (x) 

+J a »(x)A;(x) +f (x)C a {x) + C a (x)C(x)}} (4.13) 
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where 



C eff = ~F a ^F^ + ^m 2 A a ^A« - ±-(d^A*) 2 - d»C a VfC b (4.14) 

is the effective Lagrangian given in the general gauges. In the Landau gauge (a — > 0), The 
Lagrangian (4.14) goes over to the one given in equation (3.15). When the mass m tends to 
zero, equation (4.14) is immediately converted to the Lagrangian encountered in the massless 
gauge field theory. 

Now let us turn to show that the quantization described above is equivalent to the quanti- 
zation by the Faddeev-Popov approach. According to the procedure of the latter approach, 
we need to insert the following identity [14] 

A[A] J D{9 a )5[d»Al+a\ e ] = 1 (4.15) 

into the vacuum-to-vacuum transition amplitude, obtaining 

Z[0] =< 0+ I 0- >= 1 J D(A«,\ a , e a )A[A}5[d»Al + a\ e }exp(iS x ) (4.16) 

where S\ is the action given by the Lagrangian (4.1) 

Sx = J d 4 x{~F a ^F^ + l -m 2 A^A1 - l -a{\ a ) 2 ] (4.17) 

The delta- functional in equation (4.16) implies 

d^{A 6 )l + a{\ e ) a = d»Al + a\ a = (4.18) 

where {A e )" and (A e ) a represent the gauge-transformed vector potentials and Lagrange 
multipliers, (A 6 )^ = A^ + 5A^ here SA^ was denoted in equation (1.9) and (A e ) a = \ a + S\ a . 
Equation (4.18) holds naturally because the constraint condition (4.2) is required to be 
satisfied for all the field variables including the ones before and after gauge transformation. 
The SX a may be determined by the requirement that the action S\ is to be gauge-invariant 
with respect to the gauge transformations of A a ^ and A a . By this requirement, it is easy to 
find 

SSx = J d i xd^Al{8\ a - m 2 e a ) = (4.19) 

where the condition in equation (4.2) has been considered. From the above equation, noticing 
Qfj-j^a ^ o in the general gauges, we see, it must be 

S\ a = m 2 9 a . (4.20) 

Here we see that in the case of massless gauge fields, <5A a = so that equation (4.19) becomes 
a trivial identity. When the gauge transformations given in equations (1.9) and (4.20) are 
inserted into equation (4.18), we may obtain a constraint condition which is identical to that 
denoted in equation (4.6). Hence, from equation (4.15) we get A[A] = detM[A] in which 
the matrix M[A] is completely the same as given in equation (4.8). It is easy to verify that 
the determinant det M[A] , the integration measure and the action in equation (4.17) are all 
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invariant with respect to the gauge transformations of the functions A a ^ and A a . Therefore, 
when we make the gauge transformations: A a ^ — ► (A~ e )^ and A a — > (A~ e ) a to the functional 
integral in equation (4.16), the integral over 9 a {x) , as a constant, may be factored out from 
the functional integral over A® and A a and put in the normalization constant N. Thus, we 
have 

Z[0] = /" £>(A° , A a ) dot M [A]<5[c^ + a\}exp(iS x ) (4.21) 

On completing the integration over A a , we obtain 

Z[0] = 1 y D^detM^ea^iS - ^ / d 4 z(^) 2 } (4.22) 

where S is the massive Yang-Mills action. By making use of the representation of det M[A] 
as shown in equation (3.12) and introducing the external sources, we exactly recover the 
generating functional written in equations (4.13) and (4.14). 

In the end, we would like to emphasize that under the Lorentz condition which has been 
introduced into the Lagrangian in equation (4.3) and the transition amplitude in equation 
(4.16), the gauge-invariance of the action shown in equation (4.4) and equation (4.19) was 
formulated only for the infinitesimal gauge transformations denoted in equation (1.9). The 
quantized result shown in equations (4.13) and (4.14) was derived just by the use of such 
gauge transformations. In the Landau gauge (a = 0), this result is exactly the same as 
that obtained in section 3 by the quantization performed in the Hamiltonian path-integral 
formalism for which we only need to calculate the classical Poisson brackets without concern- 
ing any gauge transformation. This fact indicates that to get the correct quantized result 
by the methods formulated in this section, the infinitesimal gauge transformations are only 
necessary to be taken into account. This implies that in the physical subspace restricted by 
the Lorentz condition, only the infinitesimal gauge transformations are possible to exist. 



V. CONCLUSIONS AND COMMENTS 



The new features of this paper consist in: (1) The Lorentz-covariant quantization is suc- 
cessfully achieved in the Hamiltonian path-integral formalism; (2) The result of this quan- 
tization is confirmed by the quantization performed by the Lagrange multiplier method in 
the Lagrangian path-integral formalism. The latter method is proposed first in this paper 
and shown to be equivalent to the Faddeev-Popov approach; (3) The established quantum 
massive non-Abelian gauge field theory can be straightforwardly converted to the massless 
theory in the zero-mass limit in contrast to the previous quantization in the Hamiltonian 
formulation [1,2,5] for which the zero-mass limit does not exist. The essential points to 
achieve these results are clarified in this paper. They include: (1) The Lorentz condition, as 
a necessary constraint, must be introduced initially and imposed on the massive Yang-Mills 
Lagrangian so as to give a complete formulation of the dynamics of massive gauge fields. 
That is to say, in the whole space of the full vector potential, the massive non-Abelian gauge 
fields must be considered as constrained systems and the massive Yang-Mills Lagrangian 
itself is not complete for describing the dynamics; (2) The Lorentz condition may be in- 
corporated into the massive Yang-Mills Lagrangian by the Lagrange multiplier method so 
that each component of a vector potential acquires its canonically conjugate counterpart. 
This makes the Lorentz-covariant formulation of the Hamiltonian path-integral quantization 
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become possible; (3) In the physical subspace restricted by the Lorentz condition, the action 
other than the Lagrangian for the massive gauge field is gauge- invariant. Since the action is 
of more fundamental dynamical meaning than the Lagrangian, the gauge-invariance of the 
action implies that the dynamics of the massive gauge field is gauge-invariant. Therefore 
, the quantum massive non-Abelian gauge field theory may also be set up on the basis of 
gauge-invariance; (4) In the physical subspace, only infinitesimal gauge transformations are 
possibly allowed. This point explicated clearly in the end of the former section was already 
pointed out in reference [14]. In the reference, after introducing the identity denoted in 
equation (4.15) into the vacuum-to vacuum transition amplitude, the authors said that " We 
must know A [A] is only for the transverse fields and in this case all contributions to the last 
integral are given in the neighborhood of the unite element of the group". By this point, it 
can be understood why in the ordinary quantum gauge field theories such as the standard 
model, the BRST-transformations are all taken to be infinitesimal. 

Now we are in a position to comment on the previous works for the massive gauge field 
theories. Originally, the theory was set up from the massive Yang-Mills Lagrangian alone 
which is viewed as complete for describing the field dynamics [15-17]. To overcome the gauge- 
non-invariance of the mass term in the Lagrangian, the Stiickelberg Lagrangian in which 
the mass term is written in a gauge-invariant form by introducing the extra Stiickelberg 
field functions was proposed and widely chosen to be the starting point to establish the 
quantum theory [7-11]. All these theories were considered to be nonrenormalizable due to 
the appearance of the unphysical gauge degrees of freedom in the theories which arises from 
the finite gauge transformations. We would like to point out that according to the basic ideas 
stated before, the Lorentz condition should be , as a primiry constraint, introduced initially 
into the theories and thereby only the infinitesimal gauge transformations are necessary to 
be taken into account. Based on this point, let us analysize the previous works mentioned 
above. In reference [15], the authors proved an equivalence theorem by which they gave a 
Hamiltonian derived from the massive Yang-Mills Lagrangian by introducing an auxiliary 
Stiickelberg field functions. When making an unitary transformation to the Schrodingcr 
equation they derived, the mass term in the Hamiltonian becomes dependent on the auxiliary 
field and contains an infinite number of terms in its expansion of power series which lead to 
bad unrenormalizability. In reference [16], the equivalence theorem was given in the form 
of S-matrix. The author also introduced the Stiickelberg field and used it to make a finite 
gauge transformation to the fields involved in the theory. As a result, the mass term in the 
S-matrix contains an exponential function of the auxiliary field which gives rise to an infinite 
variety of distinct primitively divergent graphs that can not be eliminated by the introduced 
rcnormalizability conditions imposed on the gauge transformation (see equations (39) and 
(40) in the reference). It is noted here that when infinitesimal gauge transformations are, 
as they should be, considered only in the theory, it is easy to find that the gauge degrees 
of freedom are reduced and, particularly, the introduced renormalizability conditions are 
perfectly satisfied. In reference [17], the author found a relation by which any full vector 
potential may be represented as a gauge transformation of the physical transverse vector 
potential and tried to separate the gauge degrees of freedom from the transverse ones as it 
is able to be done for the renormalizable theory of the massive neutral gauge boson coupled 
to a conserved current. He eventually failed to do it because in the non-Abelian case the 
coupling between the both degrees of freedom does not vanish upon integration and thereof 
he concluded that the theory is nonrenormalizable. However, this conclusion is invalid for 
the infinitesimal gauge transformations required by the Lorentz condition because in this 
case the coupling mentioned above vanishes and, therefore, the gauge degree of freedom can 
also be separated out. As a result, the generating functional, as in the case of massive neutral 
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gauge boson theory, may also be expressed through the transverse fields like equation (3.1). 
In reference [13] , the quantization of the massive non-Abclian gauge field was implemented by 
the Faddeev-Popov operation with incorporation of the Lorentz condition in the generating 
functional. Nevertheless, the incorporation of the Lorentz condition is only for the purpose 
of improving the behavior of the massive gauge boson propagator, being not thought of 
a necessary procedure. In addition, different from the procedure stated in section 4, the 
quantization does not meet the requirement of gauge-invariance so that the ghost particle 
loses a mass term in the general gauges (a ^ 0). In particular, opposite to the procedure that 
the ghost term in the effective Lagrangian is introduced by applying the infinitesimal gauge 
transformations, the finite gauge transformations are made to the effective Lagrangian. As 
a result of these transformations, the mass term depends on the parametric functions of 
the gauge group and contains various unrenormalizable infinities. As pointed out before, 
when the Lorentz condition is incorporated in the theory, a consistent procedure is that the 
infinitesimal gauge transformations are required only. In this case, the unrenormalizable 
infinities in the mass term could not appear. Let us turn to the gauge-invariant Stuckelberg 
Lagrangian [8-10] 

L = -\F^ V F; V + \m\Al u%){A a » - (5.1) 
where are the group-valued Stuckelberg functions defined as 

j P iu _ i 

w - - -(d.u-'ur = (— — r v (5.2) 

y lUJ 

in which U — e ig ^ T and u ah = igf abc <j> c . The function u>^ for each group index a contains 
an infinite number of terms when we expand the exponential function e %UJ as a series. The 
quantum theory built by the Lagrangian (5.1) was proved to be nonrenormalizable due to 
the nonpolynomial nature of the function w". It is noted that in the Lagrangian (5.1), the 
vector field is also given a mass. Therefore, if it is considered to be physical, having three 
polarized states, we should also impose on it the Lorentz condition 

d>*L>° = (5.3) 

Substituting equation (5.2) in equation (5.3), it will be found that a necessary and sufficient 
condition to meet equation (5.3) is d^4f = which according to equation (5.2) leads to 
uo^ = 0. This implies that the Stuckelberg functions are unnecessary to be introduced. 
On the other hand, if the functions are treated as free variables, not receiving any 
constraint, there should be an integral over them in the generating functional. The integral 
is of Gaussian-type and hence easily calculated. As a result of the integration, the mass term 
of the gauge fields will completely be cancelled out from the Lagrangian. Thus, we are left 
with only a massless gauge field theory. In a word, if the quantization of the massive gauge 
fields respets the basic requirements stated in the beginning of this section as was done in 
the preceding sections, the gauge degrees of freedom will be reduced and the problem of the 
nonrenormalizability will disappear. 

From the perturbative expansion of the generating functional given in equations (4,13) and 
(4.14), it is easy to derive the bare vertices and free propagators. The vertices are the same 
as those given in the massless theory. The gauge boson propagator and the ghost particle 
one are respectively, in the momentum space, of the forms 

M^k)- *d + _____| (5 .4) 
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and 



iA ab (q) = (5.5) 

These propagators show good renormalizable behavior in the large momentum limit, indi- 
cating that the power counting argument is applicable in this case to support the renormal- 
izability of th theory. The renormalizability of the theory originates from the fact that the 
unphysical degrees of freedom contained in the massive Yang-Mills Lagrangian, i e., the lon- 
gitudinal component of the gauge field and the residual gauge degrees of freedom existing in 
the physical subspace are respectively constrained by the introduced Lorentz condition and 
the ghost equation. The both constraints respectively give rise to the gauge-fixing term and 
the ghost term in the effective Lagrangian (4.14) which just play the role of counteracting 
the unphysical degrees of freedom in the massive Yang-Mills Lagrangian. This cancellation 
also ensures the unitarity of the theory. It is easy to check that for the theory of the gauge 
field coupled to a vector current such as QCD, all the S-matrix elements given in the tree 
diagram approximation are unitary because the unphysical longitudinal part of the gauge 
boson propagator gives a vanishing contribution to the S-matrix elements. For the higher 
order perturbative S-matrix elements, it can also be proved that the unphysical intermediate 
states are cancelled out with each other so that the unitarity of the theory is still preserved 
( The proofs will be reported later) . The problem arises for the theory of the charged gauge 
boson coupled to a charged chiral current as we met in the electroweak theory [13]. For this 
kind of theory , it seems that the tree-unitarity is violated if without recourse to the Higgs 
mechanism [18,19]. We do not concern this kind of theory in this paper. But, we would like 
to note that even for the charged gauge boson theory without introducing the Higgs mecha- 
nism, the unitarity of the theory may also be preserved by means of the limiting procedure 
proposed originally in reference [20]. The limiting procedure means that the gauge boson 
propagator in equation (5.4) which is written in the a-gauge will, in the limit: a — ► oo, be 
converted to the one given in the unitary gauge 

which was originally derived in the canonical quantization by making use of the Fourier 
representation of the transverse vector potential denoted in equation (1.7). Since the a-gauge 
propagator has good renormalizable behavior, one may employ this kind of propagator to 
calculate the S-matrix and then by the limiting procedure to obtain the unitary gauge result 
so as to guarantee the unitarity of the theory. In doing this, the theory given in the a-gauge 
can be viewed as a regularization of the theory given in the unitary gauge. 
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